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Abstract. 

A thorough understanding of the energy dissipation in the dynamics of wet granular matter is essential for a continuum 
description of natural phenomena such as debris flow, and the development of various industrial applications such as the 
granulation process. The coefficient of restitution (COR), defined as the ratio between the relative rebound and impact 
velocities of a binary impact, is frequently used to characterize the amount of energy dissipation associated. We measure 
the COR by tracing a freely falling sphere bouncing on a wet surface with the liquid film thickness monitored optically. For 
fixed ratio between the film thickness and the particle size, the dependence of the COR on the impact velocity and various 
properties of the liquid film can be characterized with the Stokes number, defined as the ratio between the inertia of the 
particle and the viscosity of the liquid. Moreover, the COR for infinitely large impact velocities derived from the scaling can 
be analyzed by a model considering the energy dissipation from the inertia of the liquid film. 
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Understanding the energy dissipation associated with 
particle-particle interactions is crucial for describing the 
collective behavior of granular matter [1], i.e., large ag- 
glomerations of macroscopic particles. The coefficient 
of restitution (COR), firstly introduced by Newton as 
the ratio between relative rebound and impact veloci- 
ties [2], can be used to characterize the energy dissi- 
pation at the particle level. This number provides one 
of the basic ingredients of computer assisted modeling, 
such as molecular dynamics (MD) simulation, which 
has been developed into a powerful tool to describe the 
large scale collective behavior of granular matter in the 
past decades [3, 4]. Besides the energy dissipation from 
particle-particle interactions, the dissipation arising from 
the interstitial air or liquid has to be considered when 
coping with natural phenomena such as dune migration 
[5] or debris flow [6], as well as with various industrial 
applications such as granulation process [7, 8]. 

The experience of building sand sculptures tells us that 
the rigidity of a granular material increases as a small 
amount of a wetting liquid is added. This is largely due 
to the cohesion arising from the formation of capillary 
bridges between adjacent particles [9]. The so-called wet 
granular matter behaves dramatically different from non- 
cohesive dry granular matter while agitated, with emerg- 
ing critical behavior, such as phase transitions [10] and 
pattern formations [11], being traceable to the energy or 
force scale of a single capillary bridge. In order to gain 
insights into the dynamical behavior of wet granular mat- 



ter, it is essential to explore the COR and the associated 
energy dissipation of wet impacts. A recent investigation 
reveals that the dependence of the COR on various parti- 
cle and liquid properties can be scaled with two dimen- 
sionless numbers: the ratio between the inertia of the par- 
ticle and the viscosity of the liquid (Stokes number), and 
that between the liquid film thickness and the size of the 
particle [12]. Here, further experimental results with a 
different density ratio between the particle and the wet- 
ting liquid are presented, in order to test the scaling of 
the COR with these three dimensionless numbers. 

We measure the COR by recording a freely falling 
sphere with a highspeed camera (Photron, Fastcam Su- 
per 10K) at a frame rate up to 500 fps. Using a Hough 
transformation, the center of the sphere can be tracked 
in the images and the position of the sphere can be plot- 
ted against time (see Fig. 1). Fitting parabolic curves on 
the trajectories results in crossing points, which repre- 
sent the time when the sphere hits the ground. The height 
of the crossing points varies slightly, owing to the vari- 
ous distance of the bouncing point to the camera. The 
slopes of the two fitted parabolae at each crossing point 
yield the impact velocity v; mp and the rebound velocity 
Vreb of the sphere, respectively. Hence the normal coeffi- 
cient of restitution e n = v re b/v; m p can be determined for 
every impact. Polyethylene (PE) spheres with various di- 
ameter D and a density p p = 0.94g/cm 3 are cleaned and 
pre-wetted before use. Two types of silicone oil (M5 and 
M50 from Carl Roth) with different dynamic viscosities 
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FIGURE 1. Typical trajectory of a bouncing sphere. The 
diagram shows a PE sphere of diameter D = 7.9 mm bouncing 
on a glass plate covered with a 190 jlm thick silicone oil (M5) 
film. 



77 = 4.6 mPa • s for M5 and 48 mPa • s for M50, and sim- 
ilar densities pr lq = 0.93g/cm 3 and 0.97g/cm 3 for M5 
and M50 correspondingly are used. The bottom of the 
glass container, which has an inner area of A = 100 cm 2 
is leveled within 0.03 degrees to ensure a flat surface and 
a uniform initial liquid film thickness. A sketch of the 
set-up and a more detailed description can be found in 
Ref.[12]. 

To investigate the influence of the liquid film on the 
COR, an accurate determination of its thickness is essen- 
tial. Here, the thickness is measured optically by detect- 
ing the deflection of a laser beam (see Fig. 2 a) guided 
through the liquid film. The laser beam hits the glass 
plate with an incident angle of a and a refractive angle 
j5 = arcsin(« a i r sina/ngiass), where « a ; r and n g i ass are the 
refractive indices of air and glass, respectively. Without a 
liquid film, a part of the beam is reflected at the glass-air 
interface and then reflected at a mirror at the bottom of 
the glass plate. After a few reflections at the glass-air in- 
terface (here the number of reflections jo = 4, determined 
by the length of the mirror), the beam passes through the 
bottom of the glass plate again and is detected by a CCD 
camera. 

With a liquid film of thickness % q on top of the glass 
plate, one reflection of the beam within the liquid film 
leads to a horizontal shift of the beam by a distance 

kq = 2/zii q tany, (1) 

with the refractive angle in the liquid 7 = 
arcsin(n g i ass sin/3/«ii q ) = arcsin(« a i r sina/«ii q ). 

For the typical case that the number of reflections at 
the air-liquid interface jn q is the same as that in the 
dry glass plate (Fig. 2 a), the total shift of the laser 
beam along the horizontal plane is /ji q • /i; q with /ji q = 
jo. It is also possible, as shown in Fig. 2 b, to detect 




FIGURE 2. a) A sketch of the film thickness determination 
part of the setup. The solid lines show the beam for a dry 
container, the dashed lines illustrate the deviation with a liquid 
layer on top of the glass plate. Number of shifts by passing 
through the glass plate ; g i ass = 4, number of reflections in the 
liquid film / Uq = 4. b) ,/ glass = 4, j liq = 1. c) j glass = 3, j liq = 3. 



reflected beams with a smaller number of passes through 
the liquid film (/'ii q < jo), due to partial reflections at 
the glass-liquid interface. If the filling level is too high 
and consequently the shift of the beam is too big, an 
additional reduction of reflections in the glass plate to 
./glass < jo is also possible (Fig. 2 c). In such a case, a 
negative shift of (jo — / g i ass ) ■ / g i ass has to be considered, 
where 

^giass =2/i g i ass tan/3 (2) 

is a constant value. 

Taking all the above possibilities into account, the 
reflected beams leave the glass with possible horizontal 
shifts of Jiiq ■ /uq - (jo - j g]ass ) • Z glass , which results in 

As = [j]iq ■ kq - (jo - ./glass ) ' 'glass] COS a (3) 

at the CCD. Inserting Eq. 1 and Eq. 2 into Eq. 3, the 
liquid film thickness can then be calculated by 

Jfq + 2 Uo - ./glass )%ass tan /3 

2/ liq tany 

Note that for thin films (when jo — /'glass = 0), h\i q is 
independent on the properties of the glass plate. 

Experimentally, the camera captures a series of spots 
with each of them corresponding to a certain combina- 
tion of jiiq and /glass- The brightest spot corresponds to 
the case /n q — /glass = jo, therefore it is commonly used 
for the film thickness measurement. Fig. 3 shows a com- 
parison of the film thickness obtained from V /A, with V 
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FIGURE 3. The relation between the water film thickness 
and the shift of the laser beam As at the CCD predicted from 
Eq. 4 (solid lines) and estimated with V/A, the toal liquid 
volum over the area of the container (data points). 



the volume of liquid added, and % q given by Eq. 4. With 
a certain V, the multiple spots detected give rise to mul- 
tiple As. As V increases, the shift of each spot As follows 
a linear growth with various slopes and offsets, which 
follows the prediction based on Eq. 4 with various com- 
bination of jii q and ./glass- The systematic overestimation 
of the film thickness from V/A for all data sets is pre- 
sumably due to the volume of the liquid captured in the 
meniscus of the liquid film V me n, since the systematic de- 
viation from the estimation of Eq. 4 does not depend on 
the parameters jn q and ./glass- 

Therefore, using the variation of any spot detected by 
the camera, not only the brightest one, we can measure 
the liquid film thickness % q . This additional information 
from the spots with smaller intensity gives the opportu- 
nity to expand the range of /2i; q that can be detected, pro- 
vided that the free parameters juq and ,/giass have been 
correctly determined from the first few data points. 

Figure 4 shows that the optically obtained film thick- 
ness grows linearly with the liquid volume V, giving rise 
to a slope k = l/A and an estimation of the meniscus 
volume Vmen- Therefore, an estimation of hn q with V/A 
is also appropriate, provided that the meniscus volume 
is substracted. Nevertheless, the optical way provides a 
real-time monitoring of the film thickness. This will be 
helpful in future analysis on surface waves generated by 
the impact, and in monitoring the loss of liquid due to 
evaporation. The small deviation from the linear fit for 
the silicone oil film thickness data (inset of Fig. 4) indi- 
cates a correlation between V men and /!i; q for thin films, 
which we leave to future investigations. Figure 4 also 
shows the possibility to measure the film thickness up to 
1 cm combining the information from various spots de- 
tected. The error of lm q , on the order of w 10 \im, arises 



FIGURE 4. The film thickness obtained with As and Eq. 4 as 
a function of the filling volume of water and silicone oil (inset). 
The error bars are within the size of the symbols. Linear fits of 
the data with h\\ q = k(V — V mea ) yield the same k = 0.01 cm~ 2 
for both liquids, and V men = (2.45 ±0.10) ml for water and 
(0.84 ±0.15) ml for silicone oil. 
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FIGURE 5. Normal coefficient of restitution e n against the 
impact velocity v; mp for a PE sphere of fixed diameter D = 
7.9mm at different silicone oil M5 film thickness /!i; q . 



mainly from the fluctuations of the spot intensity. At a 
larger film thickness, the spots will become too weak to 
be detected accurately by the camera. For thin films, sil- 
icone oil instead of water is preferable, because its low 
surface tension prevents the dewetting instability. Thus it 
is used for the following experiments. 

Towards a comprehensive understanding of the depen- 
dence of the COR for wet impacts on various liquid and 
particle properties, former experiments with glass beads 
have revealed that the Stokes number and the dimension- 
less film thickness h = hn q /D can characterize the influ- 
ence from inertia and viscous damping on the COR [12]. 
The Stokes number is defined as St = p p Dvi mp /9Tj. Here, 




FIGURE 6. Normal coefficient of restitution e n against the 
inverse Stokes number St -1 for a fixed h ~ 0.04. Data with 
closed and open symbols present measurements with silicone 
oil M5 and M50 correspondingly. The inset shows the same 
data, but plotted against v; mp . 



we test the scaling with PE particles, which correspond 
to a density ratio p = Pii q /p p « 1. Without the liquid, the 
COR for the PE particle on the glass plate, depending 
weakly on Vi mp for the range of interest here, is measured 
to be 0.88 ± 0.03, averaged over 368 bouncing events. 

Figure 5 shows the COR measured with an initial 
falling height of 9.5 cm and various film thicknesses hn q . 
The number of data points for each hn q corresponds to 
the number of bouncings used for extrapolating the COR. 
The error bar arises from the statistics over 8-10 runs of 
the experiment. Except for the first bouncing data (the 
one with largest t ; i mp ), the monotonic growth of the COR 
with v'i mp could be fitted with e n = e; n f(l — v c /vi mp ), giv- 
ing rise to a limiting COR e m $ at v; mp — > °° and a critical 
impact velocity v c below which no rebound would occur. 
The decrease of the COR with the increase of hy lcl clearly 
suggests the influence of the viscous damping, since the 
energy dissipated through a viscous force depends on 
the distance that the particle travels. Different from the 
case of glass beads [12], the COR for the first rebound is 
smaller than the fitted line for all film thicknesses. This 
could be attributed to the much stronger influence from 
the mass of liquid added to the particle during the first 
bouncing for the less dense PE sphere. 

Besides the liquid film thickness, the COR depends 
also strongly on the viscosity of the wetting liquid. As 
shown in the inset of Fig. 6, the COR measured with sili- 
cone oil M50 as wetting liquid yields much lower values 
than those measured with M5. A re -plot of the COR with 
St -1 in Fig. 6 denotes a linear decay, suggesting that the 
dependence of the COR on /ziiq, viscosity rj and the im- 
pact velocity Vi mp could be characterized by the Stokes 
number for fixed h. A linear fit with e n = ej n f (1 — St c /St) 



yields a saturated COR e; n f = 0.77 ± 0.04 and a critical 
Stokes number St c = 1 1 . 64 ± 1.55. Following a former 
theory [12], we could estimate e m f = e^ry — 3p7i/2e<iry = 
0.81 ±0.03, which agrees with the measurement within 
the experimental uncertainties. This agreement suggests 
that the fact that e m f is smaller than edry is due to the iner- 
tia of the wetting liquid, which provides the dependence 
on p and h. 

To conclude, the normal coefficient of restitution for 
the impact between a spherical particle and a flat surface 
covered with a liquid film is investigated for various im- 
pact velocities, particle and liquid properties. The tech- 
nique to determine the film thickness optically and the 
calibration results are described in detail. Compared with 
the former methods [12], it provides an extended upper 
limit of the film thickness, and also allows the possibil- 
ity to estimate the volume of the meniscus. At a certain 
density and size ratio, the COR is found to decay linearly 
with St -1 , which represents the scaling with the impact 
velocity and various particle as well as liquid properties. 
The parameter e m f obtained from the fitting, correspond- 
ing to the saturated value of the COR, can be understood 
by a model considering the inertia of the liquid. A more 
detailed comparison to the model with various combina- 
tions of h and p will be a focus of our further investiga- 
tions. 

We thank Laura Meissner for performing the dry COR 
measurements. We are grateful for the support from 
Deutsche Forschungsgemeinschaft through HU 1939/2- 
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